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ABSTRACT 

AN INVESTIGATION OF THE SIDE THRUST AVAILABLE FROM 
PROPELLER SHROUDING SURFACES 



by 

SAMUEL JAMES GORDON 



Submitted to the Department of Mechanical Engineering and to the Department 
of Naval Architecture and Marine Engineering in partial fulfillment of the 
requirements for the degree of Master of Science in Mechanical Engineering 
and the Professional Degree of Naval Engineer. . 



A theory and analytical method are developed for the determination of the 
circumferential circulation distribution, over an annular airfoil of arbitrary 
included angle. The resultant circulation distribution is then used to determine 
the forces acting on the airfoil surface. 

The theory is based on the "L-Method" developed by J. Weissinger, which 
consists of replacing the lifting surface by a single concentrated ring vortex and 
an associated system of trailing vortices. The method presented here is com- 
pared to other theories in the limiting case of the complete symmetrical annulus. 
It is shown to give the same result as Weissinger found using the "L-Method. " 
This, in turn, is shown to be a very close approximation to the more accurate 
theories and to experimental results. 

A sample calculation is made, based on the method presented, of the forces 
acting on a semicircular propeller shroud, to determine its usefulness as a 
steering device in lieu of the conventional ship’s rudder. The results are com- 
pared with those for a fully movable spade type rudder with area equal to that 
of the shroud and with a high lift to drag ratio. It is found that the predicted side 
thrust from the shroud is more than twice the maximum predicted for the rudder. 
Additionally the shroud is found to give a thrust augmentation equal to approxi- 
mately 10% of the propeller thrust while the rudder absorbs a minimum of 1% 
of the propeller thrust and a maximum of 30%. 



Thesis Supervisors: Justin E. Kerwin 

Associate Professor of 
Naval Architecture 

Hal L. Moses 
Assistant Professor of 
Mechanical Engineering 
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Introduction: 

The literature contains a good deal of material on the theory of annular 
airfoils and ducted propellers. Authors such as Morgan^ and Weissinger^ 
have treated the problem of a symmetrical ring airfoil at an angle attack, 
alone, as well as in the presence of other objects such as propellers and hubs. 
Weissinger has also treated the case of a complete annular airfoil with deflected 
control surfaces. There seems to be no evidence of any previous work on in- 
complete rings. 

Since the incomplete shroud seemed, from an intuitive standpoint, to 
provide a very effective means of steering a ship, I constructed a small self- 
propelled model and fitted it with a semi-circular propeller shroud. The 
model is illustrated in Figure 1. By rotating the shroud about the propeller 
axis the side thrust developed by the shroud could be directed so as to cause 
the model to turn to the right or left. The steering device worked very well on 
the model and provided the stimulus for this investigation. 

The theory developed here, I feel, provides a means for making a reasonably 
accurate prediction of the performance of either complete or incomplete shrouds 
with a circumferential variation in circulation distribution. 
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Conclusions: 

The results of this sample calculation indicate that partial shrouds or 
complete shrouds with controllable circulation and lift distributions can be used 
to good advantage as steering devices for fluid bom vehicles such as ships and 
ground effect machines. 

The following improvements in performance over that obtained with con- 
ventional rudders are predicted: 

1. Substantial elimination of rudder drag at all turning rates. 

2. Thrust augmentation from the shroud. 

3. Improved steering capability with little or no way on. 

4. Significant decrease in turning diameter with an equivalent size 
turning device. 

The following possible shortcomings in the application of shrouds as 
steering devices require consideration and investigation: 

1. Decrease in propeller loading necessitated by shroud induced 
cavitation. 

2. Propeller induced shroud vibration. 

3. Possibility of loss of steering capability with engines stopped. 

It is hoped that further theoretical work in this area along with experi- 
mental testing will be carried out and that the apparent benefits to be gained 
from the use of propeller shrouds for steering will be employed in practice. 




SIDE VIEW OF MODEL 




VIEW SHOWING SHROUD IN POSITION FOR TURNING 

TO THE RIGHT 



FIGURE I 
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Procedure : The method utilized for determining the forces on the shroud is that 
developed by J. Weissinger and known as the "L-Method" or "Weissinger's 
Method, ” The following basic assumptions are made: 

(a) The fluid is incompressible and baratropic. 

(b) Body forces, such as gravity, are neglected. 

(c) The flow is steady. 

(d) Viscous effects are confined to a negligibly thin boundary layer 



The method consists of the replacement of the lifting surface by a single, con- 
centrated bound vortex at the 1/4 chord point and an associated system of 
trailing vortices as shown in Fig. 2. The circulation distribution, along the 
bound vortex, is determined by requiring that the velocity normal to the mean 
line at the 3/4 chord point must be zero. Since the circulation distribution varies 
along the span of the shroud, a system of trailing vortices will be shed off which 
will form a cylindrical sheet of radius equal to that of the concentrated vortex 
ring. The boundary condition is linearized by requiring that the quotient of the 
radial velocity and the free stream velocity, at the 3/4 chord point, along this 
cylindrical sheet, be equal to the slope of the section at the 3/4 point. That is: 



Once the circulation distribution is known, the radial force and axial forces 
per unit length at any section can be found by the "Kutta-Joukowski Law" which, 
in this case, states: 



which does not separate on the shroud. 



PI: 





^ ( 1 . *• ^ 3 / 4 > 



P2: 1^(9) ~ ’ pVxp(l> 9* ^^(9) 

P3: F x (0) = pjjlrp 0-* 9» 0)4-q r ^(l, 0, 
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X 



(Thickness not shown) 

FIGURE 2 

L = Shroud length. 

© = Limiting half angle of shroud. 

9 = Angular coordinate along concentrated vortex in 1/4 chord plane. 

r l/4 = R a< fr us °f concentrated vortex in 1/4 chord plane. 

P = Pivotal point at radius equal to rj /4 in 3/4 chord plane. 

$ = Angular coordinate of pivotal point in 3/4 chord plane, 

c = Displacement of mean line relative to nose tail line, nondimensionalized 

on L. 

c ' = Slope of mean line relative to nose tail line. 

\ = x/ri /4 = Nondimensionalized axial coordinate 

p = r/rj /4 = Nondimensionalized radial coordinate 

-9 = Convergence angles of nose-tail cone 

a Q = Angle of attack in plane of symmetry. 
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For a shroud which is symmetrical about 9 = 0, the resultant radial force 
R lies in the plane of symmetry and is equal to the total lift acting on the shroud. 
R can be found by integrating the components of F r (0), parallel to the plane of 
symmetry, over the span as follows: 

© 

P4: T(0) = V 1/4 J F r (0) cos 0 d0 

Similarly, the total induced drag D, acting on the shroud, is found as 
follows: 

© 

P5: D t = V 1/4 J 0 F x (0)d0 

The problem reduces to the determination of the circulation distribution 
for a given propeller geometry and loading and a chosen shroud geometry. 

The circulation distribution is represented as a Fourier cosine series. 

00 

A-23 : ^ A n cos (n0) 

n=0 

It is to be noted here that 0 is the actual polar angle measured from the plane 
of symmetry. For an incomplete shroud the condition that the circulation 
vanish at the tips is satisfied by requiring that the sum of all terms in the series 
vanish rather than that each term vanish individually. The resultant radial 
velocity induced by the shroud, at a point on the trailing vortex sheet, is deter- 
mined as a series equation in terms of a set of coefficients and integrals. 

« 6 

itl 26i qr s = ? £ I C nm W 
n=0 m=l 
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where are coefficients involving the parameters A n and I nm are geometric 
integrals. 

Since four terms in the series representation of T(q) (Eq. A -23) will 
provide a representation of more than sufficient accuracy for our purposes, we 
will restrict our consideration to n = 0, . . . , 3. We must, therefore, evaluate a 
maximum of 24 coefficients, C mn , and 24 integrals, 1^^ for each angular co- 
ordinate ($) at which we desire to determine q rs . These can then be combined, 
using Eq. A-26, to form an equation in terms of the coefficients, A n , of Eq. 
A-23. If we let n = 0, ...» n^^., ^ 3, we must evaluate q rs at (n max - 1) 

values of $ by using Eq. A-26. 

Knowing the propeller velocity field for a given propeller geometry and 

loading, and knowing the shroud geometry we can evaluate q r at the same 

s 

(n m ax “ 1) values of $ by Eq. P-1. In the case of an incomplete shroud we have 

an additional boundary condition to satisfy and that is that F^ + = 0. This gives 

us an additional equation resulting in a set of n max simultaneous equations in the 

n max unknowns Aq, . . . , A n max . For a complete shroud we must evaluate 

q r at an additional point to give us n max equations, 
s 

The solution to this set of equations provides the desired representation of 
the circulation distribution and permits the determination of the side thrust 
available from the propeller shroud as well as the drag induced by it. 

Weissinger has shown, in Ref. 1, and it is shown here in Appendix C, that 
this method yields the exact value of total lift for a complete ring airfoil in the 
limit of very high or very low aspect ratios (L/ 2 ri/ 4 ) and yields a remarkably 
good approximation for arbitrary aspect ratios falling between these limits. 

It is to be noted that this method gives the correct solution only for airfoils 
of zero thickness as it takes account only of the circulation and not of the thick- 
ness. The effects of thickness must be accounted for separately; however, their 
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effect on the total force acting on the surface is small in comparison to that of 
circulation as shown in Fig. 22 , Ref. 2. (Reproduced here as Fig. A-5. ) 
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Sample Calculations: 

Consider a propeller 



D p = 12. 0 ' 

(1 -t) = 0.80 

V s = 20 KTS 

np = 0. 625 
N = 200 RPM 
THP = SHP 



seller with the following characteristics: 

SHP = 6000 
(1-W) = 0.70 

V a = V ro = 14 KTS; = 23. 65 F. P. S. 
e rr = 1.00 




Cy = Thrust Coeff. 



T 




C T = - C^ HP ) ~ . 5Q ; p = i. 99 slugs/ft 3 
TTRp P(V/ 



2(3750) x 550 



7. 5 x 5. 5 x 10 6 



tt(6. 0) 2 (1. 99) (23. 65 ) 3 



rr(3. 6) (1.99) (2. 365 ) 3 x 10 4 



Assume the circulation distribution over the propeller is constant; radial and 
axial velocity fields for this case are tabulated in Tables 3 and 4, Ref. (8). 
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Shroud Geometry. 

Consider a shroud of constant chord, camber and convergence angle, and 
assume, for the purposes of this calculation, that the velocities induced by the 
propeller, at the shroud vortex sheet, are equivalent to those induced on a 
cylinder of radius Rp. Take: 



Referring to Table 3, Ref. (8) at ^- = 1. 0 ; -0. 4 (shroud leading edge) 




Choosing NACA a = 0. 6 meanline as tabulated in Ref. (7) 



Oj ■ 2. 58“ ; c ' I x = i- = + 0. 1827 



E I x= = °' 0737 





Referring to Table 4, Ref. (8): 




= - 0. 086 ; = -0. 086 (1. 385) v, 
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FIGURE 3 



For ideal angle of attack: 



Op = -p + = tan 



-1 




tan 



-1 



SO. 119 'N 

vttttv 



Op = - P + cc^ = 6.12° 



-p = 6. 12° - 2. 58° = 3. 54° = 0. 620 Radians 



The shroud will have the following characteristics: 



r 1/4 = 6. O' + (0. 737) (9. 6) + 0. 62 (9. 6/4) = 6. 857* 



L = 9.60' 
0 = 90° 



Mean line: N. A. C. A. a = 0. 6 
-p = 3. 54° 



, L/2 4.80 

"■3/4 rx/4 " 057 



0. 70 
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Boundary Conditions for Determining Thrust: 



(1) at 0 = ± ®;r (0) = o 

(2) at ? - !? -■ J* 1 ? = c 1 - (-P) = 0. 1827 - 

q xp + v ~ 



\s * faxp + 1207 > " Irp 

From Table 3, Ref. (8): 



— = 0. 155 

Vo>C T 



q^p = 0. 155 (1. 385) v m = 0. 215 
From Table 4, Ref. (8): 



^rp 



-0. 136 



q rp = "°* 136 385) v oo 

q = (0. 147 + 0. 188) v = 0. 335 v 

^rs co oo 



q rg = 0. 335 (23. 65) = 7. 93 F. P. S. 

Take: 

§! = 0 



0. 062 



$2 = tt/4 
§3 = tt/3 
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The set of simultaneous equations determining the four coefficients in the 
cosine series representation of the circulation distribution are: 



0. 025 A 0 + 0. 0741 Aj + 0. 0963 A 2 + 0. 1739 A 3 = 7. 93 
0. 0235 Aq + 0. 0422 A 1 + 0. 0198 A^ - 0. 0450 A3 = 7. 93 
0. 0214 A q + 0. 0248 Aj + 0. 0165 A^ - 0. 153 A3 = 7. 93 

1. 0 Aq + 0 A^ —1.0 A 2 + 0 A3 — 0 

Solving simultaneously: 

Aq = A 2 = 338.0 
Aj = -249.2 
A 3 = -84.6 

= 338 - 249. 2 cos (0) + 338 cos (20) - 84. 6 cos (30) 
By the Kutta-Joukowski Law 



'(oT) 



. _ p r r Il/2 

H COS Ctrp 1/4 J.^2 



Tg cos 0 d0 



By Ref. (8): 



^2L 



V C T 

OO 1 



X=0 



= 0. 10 ; = 0. 1385 V a 



L( aT) = -188,500 Lbs. 




is positive when directed radially outward, perpendicular to V tota j 



-17- 



Determination of the Induced Drag: 

D^(aT) = - T(0) tan 

D i (aT) = + 188,500 (0.045) = 8,480 Lbs. 

The induced drag is in a direction parallel to V tota j 




-T(0) = cos <Xp + D. (aT) sin aT 

~T(0) = 188,500 (0. 9943) + 8480 (0. 1067) 
-T(0) = 187, 500 + 940 = 188, 440 
"k(aT) = 188, 500 Lbs. 

-Dj(0) = -L^ aT j sin (aT) - D^aT) cos aT 

-Dj(0) = 20, 100 - 8, 400 = 11,700 Lbs. 
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Form Drag: 

For an N. A. C. A. 66-210 section 
By Ref. (7): 

C D © a = 0 = 0. 0094 
where 

D f 

C D | _ = 

1 P( V co + ^xp) 2 n r l/4 L 

D f = 0. 0094 (1. 0) [(1. 215) (23. 65)] 2 tt (6. 85 7) (9. 6) 

D f = 9. 4 (0. 826) tt (6. 857) (9. 6) 

D f = 1600 Lbs. 

Taking as acting in the +x direction 
Dtotal = Dj + Df =-11 , 700 +1600 

D total = ~ 9 ’ 100 lhs - 

Final Predicted Thrust from this Shroud: 

Ahead thrust = 9, 000 Lbs. 



Side thrust = 188,000 Lbs. 
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Comparison with a Conventional Rudder: 

Assuming that a/6 = 5/7 where: 
a = angle of attack 
6 - rudder angle. 

For a rudder of Aspect Ratio = 3 } for the best L/D ratio and of area equal to the 
area of the shroud: 

A = area = (6. 86 + 0. 48 ) tt (9. 6) = 221 Ft. 2 
AR = Aspect Ratio = (2S)^/A = 3 
S = Span 

S 2 = | (221) = 166 ; S = 12. 88' 

p-_k_ p _ D 

\->T O '-n w O 

(p/2)\TA (p/2)v;A 

By Figure 62, Ref. (9): 



6 


a 


6 -a 


Cl 


L 


CD_ 


D 


0 


0 


0 


0 


0 


0. 01 


930 


22.4° 


16° 


6. 4° 


0. 830 


77, 000 


0. 096 


8, 900 


28. 0° 


20° 


8. 0° 


0. 980 


91,000 


0.163 


15,120 
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T q = Side Thrust = L cos (6 - a) - D sin (6 - a) 
D q = Drag = L sin (6 - a) + D cos (6 - a) 

6 = 22.4°: 

T 0 = 76,000 
D q = 17,400 
6 = 28°: (Breakdown) 

T q = 88,000 ' 

D 0 = 27,500 
6 = 0 : 

Tq = 0 
D 0 = 930 
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SYMBOLS: 

A n 

c 

c' 

C D 



C L 



^nm 



D 

D f 

Di(a) 

E(cpp,k) 

F r 

Fx 
f»f» k 

w 

k 2 

K(k) 

l 

L 

1(a) 



Coefficient in series representation of circulation distribution 
Nondimensionalized height of mean line above the nose tail line. 
Slope of the mean line relative to the nose tail line 



Drag coefficient = — * — ? «- 

" R P A P v 2 



Lift coefficient = 



T(0) 



P y « r l/4 X 

Coefficient in series representation of induced radial velocity. 

T 

Thrust coefficient = ^ j 

" R P7 pV a 

Position vector from vortex element to pivitol point. 

Form Drag 
Induced Drag 

Incomplete elliptic integral of the first kind 
Radial force per unit arc length 
Axial force per unit arc length 
Unit vectors in Cartesian coordinates 

Integral occurring in the series representation of induced radial 
velocity. 

Modulus in elliptic functions 

Complete elliptic integral of the first kind 



X 2 + 2 



= constant for a particular pivotal point 



Chord length 
Lift force 



N 



Propeller R. P. M. 
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P 

q 

% 



^r p 

fr a 

q T 

%cp 

^xs 



—4 —4 «4 

r,cp,k 

R 

SHP 



t 



t 



THP 

T (e) 

v a 

V„ 



v l* v 2» v 3 




w 



x,§ 

°i 
- a o 



(3 

Y 



An arbitrary field point where the velocity is to be determined. 
Induced velocity. 

Velocity induced by the bound vortex distribution. 

Radial velocity induced by the propeller. 

Radial velocity induced by the shroud. 

Velocity induced by the trailing vortex distribution. 

Axial velocity induced by the propeller. 

Axial velocity induced by the shroud. 

Radial coordinate. 

Unit vectors in polar coordinates. 

Resultant force acting on the shroud. 

Shaft horsepower. 

/ 

Dummy variable = cos v 
Thrust deduction. 

Thrust horsepower. 

Side thrust. 

Speed of advance through medium. 

Free stream velocity. 

Functions used in evaluating elliptic integrals. 

“4 —4 

v ro + q = axial velocity of propeller stream. 

"r 

Wake fraction. 

Axial coordinate. 

Induced angle of attack at the 1/4 chord point. 

Angle of attack measured from incoming stream to shroud axis in 
plane of symmetry. 

Convergence angle of mean line. 

Circulation per unit length (positive in the positive coordinate 
direction.) 
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r 




Circulation (positive in the positive coordinate direction). 
Nondimensionalized axial coordinate. 



E rr 

Hp 

e,$ 



0 



x 

v 

§ 

vcpp/Q.^, k) 

P 

P = r/r 1/4 

<Pp 

§ 



Relative rotative efficiency. 

Propeller efficiency. 

Polar angle. 

Limiting half angle of shroud. 

x/r i f 4 = Nondimensionalized axial coordinate. 

Angular variable of integration = tt - (§ - 0). 

Axial coordinate. 

Incomplete elliptic integral of the third kind. 
Density of fluid. 

Nondimensionalized radial coordinate 
Amplitude of angular argument in elliptic functions. 
Angular coordinate of the pivotal point. 



Subscripts: 



B Due to bound vortex 

n, m,p Indices 

p Due to propeller 

r In radial direction 

s Due to shroud 

T Due to trailing vortices 

x In axial direction 

G,cp In tangential direction 

1/4 At 1/4 chord point 

3/4 At 3/4 chord point 
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APPENDIX A 

DEVELOPMENT OF THE INTEGRAL EQUATION FOR THE 
VELOCITY INDUCED BY THE SHROUD 




Consider the control volume enclosing the incremental arc of the bound 
vortex distribution lying between the polar angles 0 and 0 + d0 of Fig. Al. 
Helmholtz' vortex theorem states that vortex filaments cannot end in the fluid. 
This requires that the circulation around the vortex filaments leaving the 
control volume must equal that of the vortex filaments entering. This theorem 
implies conservation of vorticity and, in terms of the circulation states that: 

r(0 + A0) + Y x (0)r 1/4 A0 = r(0) 

dr fl ( 0 K 

** + Y x (e) r i/4 *0 = r < e > 



Therefore: 
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A-l: 



Y x (0) = - 



i dr(e) 

r l/4 de 



II. Velocity Induced at a Field Point by the Shroud Vortex System : 

By the law of "Biot-Savart": 

A-2 : 

<S * r < ds XD> 

4 " |D| 3 

where: dq = incremental velocity induced at a particular field point or 

"pivotal point". 

3s = directed incremental element of the vortex line. 

T = magnitude of the circulation around the vortex element 3s. 
— *■ 

D = vector displacement from the vortex element to the 
pivotal point. 



IntegratingA2, it follows that: 



A -3: 



r * )r i di xi 



and 

A-4: 



q T (r,$,X) = 4^" J H J Y x (9)r l/4^x X P) 
-0o ° 



de 



If the pivotal point coincides with a point on the vortex system, one or both of 
the integrals must be evaluated as the Cauchy principle value to exclude the 
singularity at the pivotal point. 
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ni. Integral Equations for Velocity Induced by the Bound Vortex Distribution: 




FIGURE A2 



"P" is an arbitrary pivotal point where the induced velocity is to be deter 
mined, p has cartesian coordinates (x, y, z) and polar coordinates (r, $,x). 

By the law of Biot-Savart: 



A -.4: 




4tt 



r r l/4 



3exp 
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In cartesian coordinates: 

A-5: 

de = d 0 [Of- sin 9 j*+ cos 9 k] 
r = r [ 0 "i + cos §"j + sin $ k] 

A- 6 : 

D = (x)t + [r cos $ - ri /4 cos 0 ] J + [r sin $ • rj /4 sin 0 ] k 

39 X D = -d0 £sin 0 (r sin $ ~ r 1/4 sin 0) + cos 9 (r cos $ - rj /4 cos 0 
+d 0 [x cos 0 ] j 
+d 0 [x sin 0 ] k 

| D | = [x 2 + r 2 + rJ /4 - 2r r 1/4 cos ($ -9 )] 1/2 

Nondimensionalizing: 

A- 7: 

— — = p ; -2- = X 

r l/4 r l/4 

A- 8 : 

|°| = r 1 / 4 [X 2 +p 2 + I - 2 pcos ($- 0)] 1/2 

A -9: 

3© XD = ri /4 "[” [ P cos ($-0) - l]d0 i+ \ cos 0 d0^ + X sin0 d0 k j- 

Substituting in A-7, A- 8 , A-9 into A-4: 

rln = - r (9) [p cos (i -0) - 1 3 de -f 

bX 4 lTr l/4 U2+p2 + l - 2 pcos(4-8 )] 3/2 

_ r(0 ) X cos 0 d 0 1 

4nr 1 / r 4 [X 2 +p 2 + 1 - 2 p cos ($ -e )] 3 ^ 2 
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F(8) X sin 8 d9 

4 TTr i/4 [X 2 +p 2 + l -2p cos ($ -8)] 3/ ^ 2 





d %r r - 



T(8 ) X cos 8 cos $ d8 



4n r i/4 [X^ + p 2 + 1 - 2 p cos (§ - 9 



j T(8) X sin 8 sin $ d8 

4tt r ^4 [X 2 + p2 + 1 - 2 p cos (§ - 6)]^^ 



r 




Xr 

4TTr l/4 



F(9) cos ($ - 8) d8 
[X 2 + p 2 + 1 - 2p cos ($ -9)]^^ 



Integrating over the shroud: 



A-10: 



q br 4TTr 1/4 J [x2 + p 2 + l - 2p cos ($-8)] 3/2 



Xr 



J 



© 



r(e) cos(§ - e) de 



Similarly: 

A-ll: 

‘Ikv = 



'hr l 



© 



r(8)[pcos ($-8) - l]d 8 



bx " 4 " r l/4 J .0 j>2 + p2 + l-2pcos (S-e)] 3/2 

3? 



dq b * = % 



d$ 



d %^ 



_ j -r(9 ) X cos 9 sin $ d8 

[4TTr 1/4 [X 2 + P 2 + 1 - 2p cos ($-8)] 3 ^ 2 

+r(8) X sin 8 cos cp d8 

4n r i/4[X 2 + p 2 + 1 - 2p cos ($ - 8)]^^ 



>• cp 
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A-12: 




-W r 0 r(e)sin(^-e)de 

4 rrr l/4 J -0 l\ 2 + p 2 + l - 2 p cos (7 - 0 )] S/2 



Note: q, ^ is of little interest in this problem and is included only for the sake 
of completeness. 



IV. Integral Equations for the Velocity Distribution Induced by the Trailing 
Vortex Distribution: 



Y 




— &■ 



^ — e 



ir 

z 




FIGURE A3 
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By A-l: 



Y x (0) = " 



1 

r l/4 



^9 



By the law of Biot-Savart: 





df X 5 

I d| 3 



A-13: 





df X 5 
|D| 3 



D = (x - §)i + (r cos $ ~ *i /4 cos 9) j 
+ ( r sin $ " r i /4 si n 0)k 

df = d§ 1 

df X 0 = (0)f - d§(r sin $ ~ x \/4 sin 9)7 
+ d§ (r cos $ " r i /4 cos 9)k 

| D | = [( x -§) 2 + r2 + r i /4 “ 2r r l /4 cos @ “ 0 )] ^ 

Nondimens ionaliz ing: 

By A -7: 




A -14: 

Define g 
r l/4 



e 
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A-15: 

|D| = r 1 / 4 [(X - e) 2 + p 2 + 1 - 2p cos (§ -0)] 1//2 

A-16: 

df X D = - r 2 / 4 ^de (p sin § - sin 0)7 - de (p cos $ - cos 0) k 

Substituting A-14, A-15, and A-16 into A-13: 
dq^oc = 0 

n 1 ( dr(0)^ (p sin $ - sin 0) ded0 7 

T 4TTr l/4 ^ dG ^ [(X - e) 2 + p 2 + 1 - 2p cos (§ -0)] 3/2 



_ 1 f dr(0) \ (p cos $ - cos 0)ded0 k 

^ 4nr l/4 ^ dQ ' [(x - e ) 2 + p 2 + 1 - 2p cos (§ - 0)f^ 2 



dq Tr dq T # | r | dq T 



[cos $ j + sin $ k] 



dq Tr r 4nr 



1/4 



dr(0) [+sin ($ - 0)]de d0 r 

d0 [(X-e) 2 + p 2 + 1- 2p cos (§ -0)] 3//2 



dq T * = dcL * = dq„ • [-sin $ f+cos $ k] 

l<*p| 1 



Hn - _ 1 dr(0) [P - cos ($ -0)] de d0 cp 

dq T^ - - 4nr 1/4 d0 [(x - e) + ,2 + , . 2p cos (i . 0) f/2 



integrating over the shroud: 



A-17: 



,0 to 



Tr 4nr l/4 4 J o 



IT fSSl 

-© J o d0 



sin (§ - 0) de dQ 

[(4 -e) 2 + p 2 + 1 - 2p cos ($-0)] 3//2 
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A-18: 

0 00 dr 

a = cp P P 01 (9 ) [p - cos ($ - 9)] de dQ 

T§ 4tt x i /a J J d0 ny ">2 . ~2 _ i "V r\ ~ ~ ' /i 

1/4 -q o IIX - e) -f p + 1 - 2 p cos (§ - 0)J 

A-19: 

cjfjOC - 0 

Considering A -17: 

5tt ■ J isF 1 sin <*- 0) } ® 



where k 2 Hf? + l- 2p cos ($ - 0) can be treated as constant when integrating 

o 

the integral in brackets (k ^ 0). Integrating the bracketed integral by parts: 
Let u = X - e ; du = -de 
LIMITS: e = 0 =) u = X 
e = °° =) 11 =“® 



00 




A-20: 



q 



Tr 



de _ P _ du _ u 

C(X - e) 2 + k 2 ) a/i ' J X ' [ U 2 + k 2 ] 3 ^ = 
v ’ = A. + — 2 

k 2 (X 2 + k 2 ) 1/2 k 



r A T 0 dr (9) sin ($ - 0)d0 

4nr l/4 *1-® dQ [p 2 + l - 2p cos ($-0)] [X 2 + p 2 + l - 2pcos(4-0)] 1 ^ 2 

.j r P 0 dr (0) sin (§ - 0 ) de 

4TTr l/4 ••© d9 [p 2 +1 - 2pcos (§-0)] 
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V. Complete Integral Equations for the Velocity Induced by the Shroud: 
Combining A- 10 and A-20: 

% = %r + *Tr 

A-21: 



_ © 

q = f r(e) cos ($-9) de 

r s 4 nr 1/4 J_ @ [ X 2 + P 2 + 1 - 2 p cos ($-0)] 3/2 

! \ r £® dr (9) sin (cp - 6) dQ 

4TTr l/4 [p2 + l -2pcos (§ - 6)] U 2 + p 2 + l-2pcos($-0)] 172 



, r I 0 dr (9) sin ($ - 8 ) d6 

4rr r i/4 dQ [p2 + \ - 2p cos ($ -0)] 

By A-ll: 

q_ =n - r @ r(e)[p cos (§-e) - i]de 

%x 4nr 1/4 J. @ [x 2 + p2 + i I 2p cos ($. 0 )f/ 2 
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APPENDIX B 

SOLUTION OF THE INTEGRAL EQUATIONS FOR A PIVOTAL 
POINT COINCIDENT WITH THE TRAILING VORTEX SHEET 

I. Reduction of Eg. (A - 21 ) to Inegrable Form 

For a pivotal point on the trailing vortex sheet, r = ri/4, p = 1 , and the 
integrals involving dr(g)/d0 which arise from the trailing vortices must be 
evaluated as Cauchy Principle Values. 

A -22: 

- = I* 0 r (Q) cos ($- 9 )d 9 

CiXs 4nr l /4 J-0 [X 2 + 2-2 cos (§ - 9 )] 3/2 

0 

X r L dr (8) sin ($ - 9 )d 9 

8rrr l /4 de [ 1 - cos (§ - 0 )] [X 2 +2 - 2 cos ($ - 0 y ? /2 

I r £, r (8) sin($ -0)d8 

8TTr l/4 de [1 - cos ($ - 0)] 

Assuming that the circulation distribution can be represented as a com- 
plete Fourier series: 

00 

r<e> = Y A n cos (n0) + B n sin (n0) 
n=o 

The boundary conditions on this representation are: 

a. T(0) is an even function of 0 for a shroud symmetrical about the 
0=0 plane. 

b. T(0) = r(-@) = 0 (0 7^ tt); that is the circulation goes to zero at the 
tips of the shroud. 

By (a) B n 



= 0 for all n. 
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A-23: 



r(0) = ^ cos (n0) 



n=o 



dr(8) = 
de 



■ ^ nA n sin (n0 ) 



n=o 



A-24: 



Define 



X 2 + 2 



= t 



Substituting A-23 and A-24 into A-22 and summing over the integrals: 



A -25: 



rX 



1 A„ 



n=o 



.0 



- , r w cos (n0) cos (§ - 0)d0 

rs 8/2nr 1/4 U - cos ($ - 0)] 3/2 



r X 



00 



L ^ © 

n=o £ sin (n0) sin ($ - 8) d0 



# 



8/2-nrj^ _© [1 -cos(§ -0)] [t-cos(§ -0)] 



172 



r 

- . n= j? 



00 

I nA n 



.0 



8rr r l/4 -0 






Transformation of variables: 



n0 = n$ - n(§ - 0) 

cos n0 = cos n(§ -0)cos ncp + sin n(§ -0) sin ncp 
sin n0 = cos n($ -0) sin ncp - sin n($ -0) cos ncp 
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Define: 



v = tt - ($-0) ; dv = d0 
sin (§ - 0 ) = sin v 
cos ($ - 0) = - cos v 

cos n0 = cos n(rr- v) cos n$ + sin n(n- v)sin n$ 
sin n0 = cos n(rr - v) sin n$ - sin n(n - v) cos n$ 

since 

n(n - v) = nn - nv ; 
cos (nn) = (-1) 11 
sin (nn) = (0)™ 

cos n(n- v) = (-l) n cos nv+ (0) n sin nv = (“l) n cos nv 
sin n(n - v) = (O) 11 cos nv ~ (-l) n sin nv = -(-l) n sin nv 
.'.cos n0 = (-1) 11 [cos nv cos n§ - sin nv sin n§] 
sin n0 = (-1) 11 [cos nv sin n§ + sin nv cos n§] 

Revisedlimits of integration: 
when 

0 = +0 ; v = n - (§ - 0) = n+© -M ^ 

when 

0 = -0; v=n-(§+0) = n- ©- § = v^ 

Note: Restriction: vj» V2^n $ ^ ±0 
Further defining: 
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\ sin (n§)A 

Crf = (-l) n 

8/2 n r 1/4 



ha s 



sin (ny) cos v dy 
Ki lt + cos v] 3/2 



c n3 - H>’ 



n+ l ^ n sin <n$) 



8/2 rrr 



1/4 



t - (* ^ cos (nv) sin v dy 

n3 _ J r. , -| r,~ Tl/2 

y^ [1+cos v][£+ cos J , 



_ ,_ ix n+l ^ n cos T _ |*V2 sin (nv) sinv dv 

C »4 “ <' 1) 'n4 ‘ J.. r i "j .in , ' — 



8/2 tt r 



1/4 



Vj [ 1 + COS v]['t-+COS v] ' 



C n5 = (-') 



n+l n sin (nt)A n 



8 tt r 



1/4 



'nS*! 



_ f ^2 cos (ny) sin y dv 
[1+cos y] 



V 1 



n _ /_i\n+l n cos (n§)A„ T _ sin (ny) sinv dy 

C n6=<-'T — ^6=1 — iT+ < ^r^ 



8 n r 



1/4 



V 1 



Based on the above definitions A -2 5 reduces to: 



A -26: 

J 03 6 

^s ” r Yj Yj ^ nm Wi 

n=o m=l 



Solutions to the above integrals are given in the following equations: (n<3) 
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Integral 


ESi 


m = 1: 


ioi 


A-43 




111 


A-44 




hi 


A-45 




hi 


A-46 


m = 2: 


J 02 


A-27 




h2 


A-28 




l 22 


A-29 




J 32 


A-30 


m = 3: 


!03 


A-31 




*13 


A-32 




J 23 


A-33 




h3 


A-34 



Integral Eq. 

m = 4: Iq 4 A-47 

Il4 A-48 

1 24 A-49 

134 A -50 

m = 5: Iqj A-35 

■is A-36 

125 A- 37 

135 A-38 

m = 6 ; I Q6 A-39 

I 16 A-40 

1 2 6 A " 41 

136 



A -42 
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Note: In evaluating the integrals I n ^ and which involve elliptic functions 

the following relations are used; 



k 2 = 





1 -k 2 = 



l-l 
-C+ 1 



g = 




■A = /ITT 

k z 



1 + k" 



21 

l+l 



cpp = am up = sin 



-1 



^ (£+1) (1 +COS Vp ) 

+ cos vp) 



Note : See Figure ( A-5) for quadrant selection. 



sn(up, k) = sin cpp = sn up ; cn up = cos cpp 
dnu = J 1 - k 2 sin 2 cp 



ns u = 



snu 



tn u = sc u = 



sn u 
cn u 



sd u = 



sn u 
dn u 



nc u = 



cn u 



cs u = 



cn u 
sn u 



cd u = 



cn u 
dn u 



nd u = 



dn u 



ds u = 



dn u 
sn u 



dc u = 



dn u 
cn u 



V-, (u) = — 7 T- [E(u) - k 2 sn u cd u] 
k /Z 



"V^u) = — ~tt* [2(2 -k^) E(u) - k^u - k^ snu cdu (k'^ nd^ u + 4 - 2k^)] 



3k 



V (u) = 4(2 -k 2 ) V2 ~ 3Vi - k 2 snu cdu nd^u 
3 “ 5k 71- 
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Summary of Integral Values: 

101 = ~Y { U 1 " u 2 ‘Tn* E E < U 1> - E(u 2 )]} (A-43) 

I n = - 2 - ( — Ucppk^k) - n(cp 2 * k 2 , k)l (U 2 - U 2 ) - t[Ioi3 (A-44) 

/t+1 L J /l+l 



I 21 = 2-L 2 (u 1 - u 2 ) + 41(1 - t) [V^uj)- V 1 (u 2 )] + 2(l->t) 2 [V 2 (u 1 )-v 2 (u 2 )] 

-2^[I u ]-[Ioi] (A -45) 



I 31 = 4t 3 ( Ul - u 2 ) + m 2 (l - l) [VjCui) - V^(u 2 )3 + - -t) 2 [V 2 (u 1 ) - V 2 (u 2 )] 





+4(1 - -t) 3 [V 3 (ui) - V 3 (u 2 )3 






-2-t[i 2 i3 -3 [i n 3 -2't[ioi3 


(A-46) 


!02 " 0 




(A -27) 


COS Vl 






I 12 = 2[(t+t)J + t(-L+t)"^ ] 




(A-28) 



cos V2 



I 22 = 4{^ + 4[M^i -W]} 

(-t+t)2 L 05 JJ 



COS 



cos 



V 1 






(A -2 9) 



I 32 * {x - O 2 ~ 4U + 8£ 2 " f ) (t+tf ~ 8(t 3 + 21) (t+t)"*} 



COS Vl 



cos v2 



(A -30) 



!q 3 ~ 



1 

ft-1 




yi+T - v/T 7 ! 

/t+t + ✓ t - 1 



} 



COS VI 

cos \>2 



(A-31) 
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l l3 = 2 ^ l+t 



COS 

cos y2 



- [103] 



s At N - C0S ^1 

*23 = C“T " ~T “ 4 ) JlJr t + ^3^ 

COS Vo 



^2 



I33 = yg* j^3t^ ~ 5t “ 4tt + 8t^ + 10t + —y J / t+t 



cos 



cos 



V 1 



- [ios 3 



^2 



I04 = 0 



Ij4 = 2/-t+l [ U1 -U2 + E(u 2) - E(uj)+-^y (snuj cduj. -snu2 cdu2)j 



I 24 = 3 (^- 1) ^^+1 {2u 2 - 2u x +- j^j [E(ui) - E(u 2 )] 

+ ztt 'r^r) snu i cdu i "C nd2u 2 -^r) x 

snu2 cdu 2 ]} - 2 [I14] 



134 = 8(t- l) 2 J7+T {Vi(u 2 ) - [V 2 (ui) - V 2 (u 2 )] 

+ XTT [ V 3( u 2> - V 3 (ixi> 3 } - 2(^+1) [I24] - (41+1) [I14] 

pl +COS Vj -1 

Iq 5 ^ n Ll+cos V2 -1 

I 15 = cos Vi - cos V2 " [I05] 

9 9 

I25 = cos \>i - cos V2 - 2[Ii5] “ CI05U 

135 = 3 (cos 3 vi - cos 3 V2) - 2 [I25 3 - 3 [Iis] _ 4 £l 05 ] 



(A-32) 

(A-33) 

(A-34) 

(A-47) 
(A -48) 

(A -49) 

(A-50) 

(A-35) 
(A -36) 
(A-37) 



(A -3 8) 
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I o6 = 0 (A-39) 

Il6 = sin vi " sin V 2 ~ vj + V 2 (A-40) 

12^ - 2 sin V 2 - 2 sin + 2 (sin 2vj - sin 2\%) " (± \q) + (±v^) (A-41) 

Note : Sign is negative if yp lies in the first or second quadrant and positive 

for third and fourth quadrants. 

I 36 = \ Y 2 ' sin3 Vj) - 2[I 26 ] - [I 16 ] (A-42) 
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II. SOLUTION OF THE GEOMETRIC INTEGRALS, 1^ IN EQUATION A-26 

Integrals 1^ and I^ are elliptic integrals and are most readily solved in 
terms of Jacobian Elliptic Functions. We will defer this step until last and first 
consider integrals I n 2 > I n 3 > I n 5 » and I n ^ which can be reduced to elementary 
forms. 

Using the algebraic transformation: 



cos v = t ; sin v = /l - 1 ^ 



-sin v dv = dt 




n 



0 



cos nv s g n (t) 
1 



sin nv = f n (t) 
0 



1 





2 



3 




-4(1 -t 2 )/(l-t 2 ) 




Considering 1^: 




■ C0SV 1 v (t) t dt 
n' ' 



(Note reversal 
of limit with 
sign change) 



cos V 2 (1 -t 2 ) 1//2 [£+t] 3 / 2 
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For n = 0: f n (t) = 0 

Al *02 = 0 



For n = 1: yt) = /l -t 2 

, c ° s V ! tdt 



I 



12 " + J 



By 119, Ref. (3): 



cos V 2 [f'+t - ] 



.±n/2 



3/2 



2 


k + bx>< 4± “ ,/2 


a(a+bx) (2±n)/2 


b 2 


4±n 


2 ±n 



with a = l, b = 1, ± n = -3. 



A-28: I 12 = 2 [ (l + t) 1/2 + <&(*, + 1)" 1/2 J 



COS Vi 



cos V 2 



For N = 2: f n (t) = 2t/(l -t 2 ) 

.cos V! t 2 & 



I 



%) r 



22 "' J cosv 2 U+t ] 3/2 



2 

Integrating by parts with u = t , du = 2t dt 



dv = 



dt 



a+ty 



3/2 



V' = 



-2 



(t + t) 



1/2 



*22 " 2 



o 


COS Vi COS Vi 


2t^ 


+ 4 r t & 




cos V2 (l+t) 1 / 2 


- 


cos V 2 



By 119, Ref. (3): with a = l, b = 1, ±n = -1 
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Al29l I 22 = 4 



(£+t) 



1/2 



+ 4 



(l+t) 



3/2 



- U,l+t) 



1/2 



cos v \ 
cos \>2 



For n = 3: f n (t) = (4t 2 - 1) /l -t 2 



COS Vl 



I 32" 4 J, 



t dt 



COS VI o < 

r t^dt _ r 

COS V 2 [t+t] 3/2 ** cos V 2 U+t] 3/2 



q 2 

Integrating the first by parts with: u = t , du = 3t dt 



dv = 



dt 



(1 + t) 



372 



v = 



(£+t) 



172" 



The first integral equals 



2t c 



(l + t) 



172 



C ° S Vl + 6 J C ° S V ‘ ' 2dt 

cos v 2 



"cos V2 (-t+t) 1 / 2 



By 112, Ref. (3): 



I 



x 2 dx _ 2 (8a 2 - 4 abx + 3b 2 x 2 



7 a+bx 



15 b" 



/a + bx 



with a = l, b = 1: 



st- 



172 + if ^ 8/£/2 " 4/Lt + 3t;2 ) (l+t) 1 ^ 2 



(l + t) 

The second integral is just A -2 8 combining: 



cos 



V 1 



COS V-; 



A -30:- 



! 32 = (3I 2 -44t + 8t 2 -|) (l+t) l/2 - (8t 3 + 2^)(t + t)' 1/2 | 



COS Vj 

COS V 2 
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Consi dering I^: 



* ■ i; 



^ cos nvsinv dv 



[1 + cos v] [t + cos vf/ 2 



COS Vi 






COS V 2 [1+t] [<t+ ] 1//2 



For n = 0: gu(t) = 1 



T = f C0SVl dt 

03 ''cos V 2 [1+t] U + t] l/2 



By 56, Ref. (3): 



I 



dx 



v/u 



1 . b'/U - /KB' 

TIT b'/U + ZEE 7 



with 



A-31: 



v = 1 + 1 = a ' + b 7 t ; a 7 = lb 7 = 1 
u = £ + t=a + bt; a = t, b = 1 
k = ab 7 - a 7 b = (£ - 1) 






< tn 



' fu+i - 



t + /t-T 



COS VI 

cos 2 



For n = 1: g n (t) = t 



cos v, 

■J 



l 13 



t dt 



cosv 2 [l+t][£+t] 



T72T 
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Using v and u as defined in IQ3, t = v - 1 



l 13 



= 1 



dt 



C0S V 1 , COS vl cos VI 

(v - 1 )dt _ J* dt _ J 1 

cos V2 v/TT J cos V2 /u - J cos V2 v/TT 



= {2/l 



+ t 



/T 7 "! 



In 



fl+t - Jl- 1 



/£+t + /-l- 1 



COS VI 

cos V2 



A -32: 



I 13 = 2/ -t+t 



cos vi 



cos 



v 2 



- [^3] 



For n = 2: g n (t) = 2t 2 - 1 = 2(v - 1)^ 

g n (t) = 2v 2 - 4v + 1 



.cos v 



1 dt 



.COS Vl , .COS V 1 j, r.y 

123 = 2 J 1 - 4 J — + J 

J cos V2 /"u" J cos V 2 /V J COS V 2 vAT 



since u=£ + t= l + t-(l - t) = v - (1 - -t) 
dt = dv 



the first integral = 2 

*1 

By 111, Ref. (3): 



1 + COS Vl 7 

1 v dv 



+ cos V2 l)+v 



-Eg-. - 2(2a : bx) /ITbx 

J a+bx 3b 



with: x = v, a = (-C, - 1), b = 1 



the first integral = - (2-C, - 2 - v) 1) + v 



1+cos vl 
1 + cos V2 
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reverting back to t: 
4 



=) - J(2t-3-t) /Z+t 



COS Vj 

cos V2 



the second integral = -8 / Z + 1 



the third integral, from A^q = 



COS Vj 

cos \>2 
1 



/£- 1 



•tn 



J t+t - / - 1 



/ -t+t + / -i/ - 1 



COS Vj 
cos \>2 



Combining; 



*23 " - 4^/1 



+ t + — - — 



St+t - /£- 1 



-t - 1 ^ _A/+ 1 + / 1_ 



COS VI 



cos 



V 2 



A -33; 



*23 - (it " T ■ 4 J 



cos vi 
cos y2 



[■os] 



For n = 3; g n ( t) = 4t 3 - 3t 

= 4(v - l) 3 - 3(v - 1) 
g n (t) = 4[v 3 - 3v 2 + 3v-1]-3v + 3 
= 4v 3 - 12v 2 + 9v - 1 

By the procedure followed in solving I23'. 

1 + cosv! y 2 dv 



J 33 " 4 



a+cos V2 '/(£-!) + 



12 






1 + cos 



VI 



v dv 



+ cos V2 /(£- 1) + v 



+ 9 



f cos V 1 _dt_ _ p COS V 1 dt 
XOS V9 sz J cos y2 v'/u 
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By 112, Ref. (3): 



x 2 dx _ 2 (8a 2 
/a + bx 



4 abx + 3b^x 2 ) 

Ti? /a+bx 



with x = v, a = f, - 1, b = 1, the first integral yields 




1+cos vi 
1 + cos \>2 



reverting back to t: 

COS Vj 

cos V 2 



=) yg- |^8f. 2 - 16-t + 8 - 4f, + 4 - 4ft + 4t + 3 + 3t^ + 6tJ / -t+t 




+ lOt - 4 ft + 8-t 2 - 20f,+ 15 




COS VJ 

cos V 2 



The second integral yields, from the development of A-33: 

[COS Vi 

8(2f. - 3 - t) /ZTt 

cos \>2 



The third integral yields: 



18/Z+t 



COS \>l 

cos V 2 



The fourth integral, from A- 30 yields: 



1 p A+f “ / £ “ 1 

— to — — — 

/ £ - 1 \TU£ t" + /£ - 1 



cos vj 
cos \>2 



combining: 
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a 33 



= < 



4 [ 3t2 



+ lOt - 15t - 4 fc + 8l 2 - 201 + 30-t + 15 - 45 + 

4 



3 



SI + 1 - - ---- - 

/-t - 1 



/ £+t ~ - 1 


\ 


jf <L + 1 + /£ - 1 


1 



COS 

cos \>2 



j| .[3t 2 - 5t - 4 + 8-t 2 + 10-t + ^ 



Si + 1 - 



Si - 1 



•in 



7i+t - Jl-1 
Sl+t + / 1 - 1 



COS Vi 

cos \>> 



A -34: 



r 33 = IT ( 3t2 " 5t - 4 it + 8 1 2 + 10 1 + ^)/i+t 



COS V 



cos 



1 

v 2 



[lQ3 3 



Considering I n 5 : 



Vo COS Vi . . 

. I* cos nv sin y dy _ (* dt 

n5 - J [1 + cos v] J_ 



5 n' 



Vi 



cos 



V2 



[1+t] 



For n = 0: g n (t) = 1 

.cos V! dt 
cos V2 [1 + t] 



*05 - I 



= in(l+t) 



cos 



V 1 



A -35: 



I 



05 



In 



1+cos vi 
1 + cos \>2 



cos V 2 
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For n = 1: g n (t) = t 

I = f COS V1 tdt 
^ cos \>2 [1 +t] 



By 32 Ref. (3) 

I FFEJ = [a+bx-a ta(a+bx)] 



with x = t, a = b = 1 




[l +t - In (1+t) 



COS Vj 

cos V2 



IfS = cos vj - cos V2 - In 



1+cos vi 
1 + cos V2 



A-36: 



Il5 = cos vj - cos V2 - [l05l 



For n = 2: g n (t) = 2t 2 - 1 



*25 = 2 



COS Vi cos vi 

1 t 2 dt _ [* dt 

cos V2 t 1+t 3 cos ^ Cl+t] 



By 36, Ref. (3): 



J a+Ex = “3 [l ( a +bx) 2 - 2a (a+bx) + a 2 -Cn (a+bx)] 



with x = t, a = b = 1, the first integral of I25 =) 



2 [~ ^" " 2(!+t)+ to (1+t)] 



COS v^ 

cos V2 
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combining with -Iq5 




t-l + itaU+o} 



COS Vl 

cos V2 



I25 = (jcos^ " cos^ V2) - 2 (] 



COS \q - cos 






Li+cos v 



A-37: 

O 9 

I25 = cos^ vi - cos V2 - 2 [I15] - [I05] 

For n = 3: g n (t) = 4t 3 - 3t 

COS Vi 

I = 4 f C0S V1 t3 dt _ 3 P t dt 

33 cos V2 0+t3 cos \>2 1-1 +t] 

By the recurrence relation. Equation 2-1112, Ref. (4): 

P x n dx _ x n a x n "l dx 

J (a + bx) n^ b (a + bx) 



with x = t, n = 3, a = b = 1 the first integral of I53 yields: 
t 3 _ rt 2 dx 

3 " J jr+ty 



From I2J and 1^5: 




- 2(l+t) + to(l + t)]} 



COS Vj 

cos V2 



COS V2 

- 3 [1 + t - to(l + t)] 



cos \>2 
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combining: 



I35 = | [(cos - (cos V2>3 J -2 [cos 2 vj - cos 2 V2 



+ 



(cos vj) - (cos 



+ cos vi - 
+ cos \>2 - 



A- 38 : 

I35 = 3 [(cos V 2) 3 - 2 [i 25 ] - 3 [li 5 ] - 4 [lo 5 ] 
Considering 



v 2 



x n6 



•r 

VI 



sin (ny) sin v dy 
UT+ cos vj 



I 

-I 



COS 



cos 



V 1 f n (t) dt 

* [1+ 'l 



For n = 0 : f n (t) = 0 



A- 39 : 

106 = 0 

For n = 1: f n (t) = A -t 2 



l 16 



-r 



v i /TV 



cos \>2 



[ 1+t] 



dt 



»COS Vl 

J cos v2 




dt 



By 60 , Ref. ( 3 ): 



= E^ - 



_k_ 

2b 



~ 2 

KTEb' 



tan 



-1 




with v = 1-t = a 7 +b 7 t; a 7 = 1, b 7 = -1 
u = 1 + t = a + bt; a = 1, b = 1 
k = ab 7 - a 7 b = -1-1 = -2 
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I 



COS vl 



cos V2 



AZ 

/i + 



L-[a- 



t 2 + 2 tan' 




■,COS Vl 



1 -tJ 



cos V 2 



By the identity 



cos 






+ cos X 

- COS X 



= tan 2 (tt - x) 



-- 1 * K" - 'h) 



A-40: 



I 16 = sin vi - sin v 2 " Vj + V 2 
For n = 2: f n (t) = 2t /l - 1 2 



l 26 



cos v r cos vi (1 v 

= 2 J wr+f dt = 2 J -f4 dt 

cos V 2 cos V 2 / 1 t 



by 161, Ref. (3) 



I- 



xdx 



/Srrp 

by 171, Ref. (3): 
x 2 dx 



77^7 



J- 



/a^ - x 2 
with x = t, a = 1 



x /— 9 , a . -l x 

2 ^'X + -T s in ^ 



t-6 = 2 {0 1 +£>^ ' 2 sm- 1 COSV1 

cos V2 
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J 26 = 2(sin " sin V 2 ) + cos sin - cos V 2 sin V 2 




A-41: 



^ = 2 sin V 2 - 2 sin vj + 2 (sin 2 v^ - sin 2 V 2 ) " (±Vj) + (± V 2 ) 

Note: Sign is negative for Vp in first or second quadrants and positive for third 
and fourth quadrants. 

For n = 3: f n (t) = 4t 2 J\ -t^ - /l - t2 
Following the method used in 126* 



l 36 



= -4 



COS 

-» 




COS vi 


1 


t 2 (l -t) dt 


- r /h 


cos V2 


fl-t 1 


cos V2 /ITT 


pCOS VI 


t 3 dt 


, cos vj j.2^ _ |,cos vi /]"T£" 


J cos V 2 


/TT2 . 


COS V 2 COS V£ /1+t 



By 175, Ref. (3): 

with a = 1, x = t, the first integral yields: 



- (sin^ vj - sin^ V 2 ) + 4 (sin \q - sin \^) 
The second integral, from yields 

- 2 [l26J - 4 ( sin Vi - sin V2) 
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The third integral 




] 



A-42: 

*36 = 3 { sin3 v 2 ' sin3 v l} ' 2 [l26] - [h6] 

Considering 1^ and I n 4 : 

These integrals must be solved by use of Elliptic integrals and Jacobian 
Elliptic Functions. 

The solutions for these integrals can be found by transforming them into 
algebraic form as was done for the integrals Ij^, In3> I n 5 » and I n 5 * Upon 
transformation it is found that they all are of the form solved in series 235 of 
Ref. (5). In utilizing these solutions, the parameters will have the following 
values: 

# s 

a = 1, b = -1, c = -L, yp = cos Vp 

k 2 = a^b. = _2_. k ' 2 - i _ k 2 = ill 

K a -c l+l ’ K 1 * l+l 



S = 



/a^c* Ji+l 



; — s- = J~U+ 1 ; 1 + k 



<2 _ 21 
1+ 1 



cPp 



= am 



/(a-c)(y p -b) /(-t + l)(l + cos v p 

u p = sin " V ^“b)(v~-c r = sin ’ v — 



Up = F(cpp,k) = Incomplete elliptic integral of the first kind. 

E(cpp,k) = E(up) = Incomplete elliptic integral of the second kind. 
TT(cpp,a 2 ,k) = rr(up,a 2 ) = Incomplete elliptic integral of the third kind. 
VJ = TT - (§+ 0) ; V2 = tt - (cp - 0) 
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Note: Due to the restriction on this group of integrals that 
(a s y > b > c) yp > b = -1 

vi» V 2 7 tt 

Considering I n j: 



! nl = J 



V2 cos nv cos v dv 



vi [£ + cos v] 



3/2 



= J 



cos V 1 g n (t) t dt 



cos V 2 (t+t) A -t 2 v/X+t 



*nl = 



n COS Vl 



g n (t) t dt 



” cos V2 (-t+t)/i -t /l+t yr+t 



For n = 0: g n (t) = 1 



,.cos 

Joi = j 

cos • 



^2 

V2 



dt 



L ^ T ^ J /r^t A+t yz+i 



Considering the bracketed term: 



t = (t+t) " fe [ZTtT = 1 'n«T 



l 01 



= 1 



COS VI 



dt 



cos v 2 /I - 1 /t - (-1) /t - (-1) 



-l f 



COS VI 



dt 



J cos V 2 [t -(-£)] /I - 1 /t - (-1) /t - (-£) 



By 235, 00: 



u 

I 



dt 



^(a-g)(t-b)(t-c) 



= gF(cpp,k) 



with a = 1, b = -1, c = -l, the first integral in Iqj =) 
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/ -L+ 
By 235.01: 



- [f(cp X , k) - F(cp 2 » k)] = ~== [u x - u 2 ] 



I 






dt 



(b?c) E ^ ,k ^ 



o (t-c) /(a -t)(t -b)(t -c) 
with a = 1, b = -1, c = -l the second integral in Iqi => 
[E(cpi,k) - E(cp2>k)] 



- 21 



(l-l)JT+T 



A-43: 



101 = 7zvr { U1 " u 2 " rn [ E(u i> " E < u 2>] } 



/Z+T 
For n = 1: g n (t) = 1 



I 



11 



■r 



pCOS VI 


t 2 1 


dt 


J cos V2 


q+t) 


/TT/T+t /i+Z 



Considering the bracketed term: 

t 2 = [q+t )-^] 2 = q+t ) 2 - 2-iq+t)-H , 2 
q+t) q+t) q+t) 



= l + t - 21 + 



£ + 1 



q+t) 

= 1 ~ l [} "q + 1) 

comparing with Iq^: 



= t-l+- 



Z+ t 



] 



'll 



= 1 



cos 



V 1 



t dt 



COS V2 yi -t /t- (-1) /t - (-1) 



toi] 
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By 235. 16: 







dt 



/(a-t) (t -b) (t - c) 



gk m J 



“p 



0 



ck^ 

b 



\m 



sn-^ u j 



(1 - k 2 sn 2 u) m 



du 



with m = 1, a = 1, b = -1, c = -l 

gb m = = ~ 2 . . 

Sl+1 + 1 

The integral in I-q becomes 



-2 T U 1 (1 --tk 2 sn 2 u) du +2 r u 2 (1 - £k 2 sn 2 u) du 

/£+■! o (1 -k 2 sn 2 u) yi+l o (1 -k 2 sn 2 u) 

By 340. 01 with a 2 = k 2 , a 2 = -tk 2 



f l -=. ^ i sn iu dn = J_ 
^ 1 - k 2 sn 2 u k 2 



£(k 2 - ^k 2 )rr(cp, k 2 , k) + £k 2 uj 



= [(1 -'On(cp, k 2 , k) + tu] 



A-44: 



111 = ^ rTr(§i,k 2 ,k) - n($2» k 2 , k)~i ■ — 2 — (uj -u 2 ) - l 

/TTi L J n ~ 1 



fi + T 



Where, by 111. 06, Ref. (5): 

^E(cp, k) - (k 2 sincp coscp//l -k 2 sin 2 cpj 



tt (cp, k 2 , k) = 



k" 



For n = 2: g n = 2t 2 - 1 



* COS \>) 

I 

COS 




dt 

/(l-t)[t-(-l)][t -(-*,)] 



- 61 - 



Considering the bracketed term: 

2t 3 - 1 _ 2[(-t+t) - if - [(t + t) - i] 

[t+t] l + t 

s 2(£+t) 3 - 6jt+t}h+ 6(-t+tH 2 - 2-t 3 - Qt+t) + 1 

(£ + t) 

= 2(t+t) 2 - 6«t+t) + 6l 2 - 1 - 
= 2 1 2 + 4{t + 2t 2 - 6i 2 - 6-tt + 6i 2 - 1 - [jjqrg-j 
= 2t 2 -2tt + 2t 2 -l-[ 2 ^] 

- 2t2 - 2 K t - l+ v^)) + xdt> - 1 

Comparing with 1^ and Iq^: 

l 21 = 2 J y — — — = = - 2-tfln] - [l 0 i] 

J '/(l - 1) [t-(-l)] [t -(-£)] L J L J 

Again citing 235. 16, with m = 2 

2 j* t 2 dt 2 [ * (1 - £k 2 sn 2 u) 2 du 

J A-t) Ct-(-l)] [t -(-*)"] J o (1 - k 2 sn 2 u) 2 

_2 f (1 ~ tk 2 sn 2 u) 2 du 
o (l-k 2 sn 2 u) 2; 

By 340. 02: with a 2 = k 2 , a^ 2 = £k 2 
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f (I ■ ~^ 2 -- s - ” ) 2 - u - ^ — =-i- I" i 2 k 4 u + 2-Uc 2 (k 2 - -tk 2 )Vi + (k 2 - ^Jk 2 ) 2 V 9 1 
J n -v2 cn 2 n \2 k 4 L- 1 z J 



c 2 sn 2 u) 2 
(1 - k 2 sn 2 u) 

By 315: (since a 2 = k 2 , Vj = ^j) 



Vl(u) = [E(u) - k 2 sn u cd u] 

k /z 



V 9 (u) = — [2(2 -k 2 ) E(u) - k /2 u] -k 2 snu cdu (k /2 nd 2 u+4 -2k 2 )] 
3k /4 



V 3 (u) = 



4(2 -k 2 ) V 2 - 3Vj -k 2 snu cdund 4 u 



5k 



/ 2 



Note: Vp V 2 , are used to evaluate both ^p Igp and They are most 
readily applied by evaluating them and substituting in their values. 

f — — 2 ^ 2 - = [t 2 u + 2^(l-t)V 1 (u)+(l-t) 2 V 2 (u)j 
J (1 - -tk sn u) L 



Combining 



A-45: 



I 21 = 2t 2 (uj - U 2 )+ mi-l) [Vjfuj) - Vjfuj)] 

+ 2(1 -l) 2 [v 2 (up - V 2 (u 2 )] - 2t[l n ] - [l m ] 
For n = 3: g n (t) - 4t^ - 3t 



! 31 = I 



dt 



cos V 1 p 4t 4 - 3t 2 f 

cos V 2 *- (t+t) -* /(1-t) [t- (-1)1 [t -(-£)] 
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Considering the bracketed term: 



4t 4 - 3t 2 = t 2 (4t 2 - 3) 
('L-i-t) " (i+t) 



t 2 = [(t + t) - l\ 2 = (l + t) 2 - 21(1 + t) + l 2 



2., .72 






= 4t 3 - 4-Ct 2 + (4I. 2 - 3) [1+ t) - 21] + 



4t 3 - 4tt 2 + 4t 2 t - 3t - 4 4 3 + 3£ + ^~y^- 



= 4t 3 - 24 {24 2 - 2tt + 2t 2 - 1 - + <--3t ^ 



2 

(£ + t) 



- 4t 3 -24{24 2 -2te + 2t 2 -l- 2 ^]} + 3(4-t^)-24[l- ri - 



comparing with and I 



10 



.cos v, 

= 4 I . 



31 



t 3 dt 



J cos /(I -t) [t - <-1)1 [t - (-£)] 
considering the first integral: 



24 [^l] ' 3 [ I ll] ‘ 2 ' t L I 01 J 



By 235. 16, with m = 3: 



, f * (1 - -tk 2 sn 2 u) 3 du _ , (* 2 (1 - -tk 2 sn 2 u) 3 du 

L ~~Z 72 2 73 J ~ 72 T~3 

o (1-k sn u) o (1-k sn u) 
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By 340. 03, with a 2 = k 2 , a^ 2 = -tk 2 : 



f du a -4 I” l\ 6 u + 3t 2 k 4 (k^tk 2 ) Vj + 3 f,k 2 (k 2 - £k 2 ) 2 V, 

J (l-k 2 sn 2 u ) 3 k 6 L 1 

+ (k 2 - -tk 2 ) 3 V3J 



r a . - . ^ 2 ^ 2 . ^) 3 du = rt 3 u+3-t 2 (i--t)v 1 + 3.t(i-<t) 2 v 2 + (i--t) 3 Vo~| 
J (l-k 2 sn 2 u) 3 L J 



A-46: 



I 31 = 4-t 3 (u x -u 2 )+ 12 1 2 (l-l) [v^uj) - Vj^)] 
+ 12 ^(l-t) 2 [v 2 (u 1 )- V 2 (u 2 )] 



+ 4(1- -t ) 3 [V 3 ( U1 ) - V 3 (u 2 )j - 2^[l 21 ]-3[l u ] -2^[l 01 ] 



Considering ^ 4 : 



I 



v 2 

-I. 



sin nv sin v dv 



n 4 1 1 

V 1 [1 + COS v] [<t + COS v]2 



= i 



cos vi f n (t) dt 



cos \>2 [l+t] [-t+t ] 2 



Series 235, Ref. (5) will again be used to solve these integrals. 

For n = 0: f n (t) = 0 

A -47: I 04 = 0 

For n = 1: f n (t) = / (1 "t 2 

p cos Vj r / 1-t 2 -| dt 
14 ios vo ( 1+t )-* [£ + t] 2_ 
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considering the bracketed term: 



/T-? , /T^t /l+t _ /TT 

• L+t /l+t 



'u'l 



cos V 1 /I" 1 ? dt 

cos V 2 ^(t-(-l)] [t -(-£)] 



By 235. 07: 



■^ (t-bHV-c T dt = < a -b)«I 0 Pcd2u 



du 



with a = 1, b = -1, c = -l 



*14 



4 r u i 



r u i o -4 r U2 9 

cd z u du ■ • - J cd^ u 



fZ+1 J o 



du 



l+l o 



By 320. 02: 



I 



cd^ u du = 



— ?y [u - E(u) + sn u cd u] 
k z 



A-48: 



Ij 4 = 2/T+T[ui -U 2 + E(u 2 ) - E(uj) + ^-j-(snuicdu 1 - snu 2 cdu 2 



For n = 2: f n (t) = 2t /(I - 1^ 
cos vi 



i = r " r 2t/i -t^ i dt 

^COS V2 ^ (* +t ^ ■* [^+t]2 
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Considering the bracketed term: 



2t/l -t 2 2t/3Tt 
“^" = Sl+t 



= 2(l+t) 




- 2 



{Tt 

/T+t 




Comparing with 1^: 



^24 



COS VI 

= 2 J yT tI£r^ * - 2 [i i4 ] 

cos V2 / -t+t 



By 235. 09: 

,yp 



°b 



(a - t)(t-b) 
(t-c) 



U P 

= (a - b) (b - c) gk 2 J sd 2 u cd^ u du 



With a = 1, b = -1, c = -t, the integral in =) 

. U 1 



i 2 < 1 - J > C y^)(ra)] [ I 0 sd2 u cd2 u du ' £ 2 sd2 u cd2 u du 



By 361.27: 



sd 2 u cd 2 u du = 1 f -2k' 2 u + (l+k /2 ) E(u) 

3k 4 k /2 L 

+ k 2 (k /2 nd 2 u - 1 - k' 2 ) sn u cd uj 



since: 

1 _ (t+1) 2 1+k' 2 2-t k 2 _ 2 

k 4 " 4 ’ k /2 l-l ' k /2 (FTJ 
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'sd 2 ucd 2 udu = i!^i[-2u + Ji r E(u) + jlj- 



(jid 2 u - sn u c< 3 u] 



A-49: 



124 = 3 (t- 1) AAA {2u 2 -2uj +^x[E(u 1 ) - E^)] 

2 r, ,2 2-t \ , 

+ l+r L (nd u i "nrJ snu i cdu i 

'C nd2u 2 - Tn:) snu 2 cdu 2]} - 2 [ i h] 



For n = 3: yt) = (4t 2 - 1) A -t 2 



l 34 



COS v 



1 j~ (4t 2 - l)/l -t 2 



cos 



V 2 



<1 + t) J a+t )2 



] 



dt 



Considering the bracketed term: 



(4t 2 - 1) A - 12 = 2t (2t/l+t2) _ A -t 2 



TT+TT 



(i+t) TT+tT 



Substituting 2t = 2(1 +t) - 2 =) 



a+FTS AtATT 2 /l -t 2 

- (1 + t) (l + t ) 



Comparing with I41 and I42* 



. 4 j<=~ dt . 2 r -j . r J 

014 J cos v 2 A~fT L Z4J L J 



with t = (t + t) - t, the integral in I43 becomes 
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COS Vj 



j J( 1-t) [t-(-l)][t-(^)]dt - 

COS V9 V 



COS V 



’i Ai-t)[t-(-m dt 
cos V2 /t - (-t) 



The second integral is the same one as was solved in I 43 



.41 [ C ° S Vl Ai-t) Lt-(-l)] 
cos \>2 ' /t - (-•£,) 




Considering the first integral: 
By 235. 14: 



J (a -t)(t -b)(t - c) dt = g(a -b)(b - c) 2 k 2 J sd 2 ucd- 2 und 2 udu 



with: a = 1 , b = - 1 , c = -l, the first integral =) 

4 ^. — J (2) (t - l) 2 J sd^ucd 2 und 2 udu - J sd^ u cd^ und 2 udu 

By 361. 18: 

J sd^ ucd 2 und 2 udu = — ^ -I2 + (1 + k 7 ^ ) I4 - k '2 I^J 

k 

where I2 = Vi 

I4 = V 2 as defined in I21 
l6 = V 3 

A-50: 

I 34 = 8a-l) 2 /l+T{v 1 (u 2 )-V 1 (u 1 )+-^ T [v 2 (ui) - V 2 (u 2 )]+ ( j ^ j ) 

[v 3 (u 2 ) - V 3 ( Ul )] } - 2(1+1) [l 24 ] - (4^+1) [l 14 ] 
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APPENDIX C 

COMPARISON OF RESULTS WITH KNOWN RESULT FOR A 
COMPLETE SHROUD 

By the Kutta Joukowski Law: 

(Equation P-2) F r (9) = -p (1,9, 0) 

Taking the total lift in the 0 = 0 direction: 



TT 




In the absence of a propeller 





By A-23: 



00 




n=o 




Because of orthogonality: 
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A-51: 



T(0) = -p V ra r : / 4 Ajn 

To find the side thrust on a complete ring, in the absence of a propeller, 
we need to evaluate Al. 

For a complete shroud, of constant chord, at an angle of attack, the side 
thrust is independent of the section shape (Ref. 6). We will consider the case 
of a complete shroud consisting of a frustrum of a cone as shown in Fig. A4. 




FIGURE A4 

The shroud has a chord of L, a radius at the 1/4 chord point of rjy 4 , a con- 
vergence angle g and is set at an angle of attack, in the 0=0 plane of a(0) = a Q . 

C' (§,\) = 0 



cx(§) - -g + CCq cos $ 
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By Pi: 

A-52: 

q rs | = r V m [+B - a 0 cos $] 

l(l,cp, X 3/4 ) 

By A-26: 

co 6 

q r s * X Z C nm V 
n=o m=l 

In evaluating the integrals I nm we must abide by the restriction that 
v P V 2 ^ rr 
where 

vi = n - ($ + 0) 

V2 = tt - ($ - ©) 

Since 0 = tt, for a complete ring, 

§ 4 TT 

Vj = - $ 

\>2 - 2tt - $ ; \>2 ~ vj = 2 tt 

Note : In evaluating the elliptic functions E(cpp, k) and F(cpp,k), the following 
relations are used: 

E(-cp,k) = -E(cp, k) 

F(-cp, k) = -F(cp,k) 

E(mrr±cp,k) = 2m E(k) ± E(cp,k) 

F(mrr±cp, k) = 2m K(k) ± F(cp,k) 
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where: 

E(k) = Complete elliptic integral of the second kind. 

F(k) = Complete elliptic integral of the first kind. 

Determination of cpi and cp 2 from vj and 

/ (£ + 1) (1 + cos Vp) 

'PP = S “' 7 2(Uc« vp) 

2 sin^ $p (t+cos vp) = (£+1) (1 + cos vp) 

= 1 + t cos vp + 1 + cos vp 

(-i+cos vp) (2 sin^ cpp - 1) = 1 + cos vp 

1 + l cos vp 
• C0S<2 ‘PP ) J T +- cos Vp- 

when vp changes by 2tt, 1 2 cpp | changes by 2tt. To find the sign, assume 

1 = 2 . 

(1) let vp 1 = 0 

-cos 2cp Pj = 2 = 1 5 cos 2cp Pi = -1 ; cp Pi = -J 

(2) let vp = tt/3 

- cos 2cpp = — — = 0. 8 

2 + i 

2cpp 2 = tt - 37° cpp 2 cpp 2 = j ~ 18° 

(3) let vpg = rr/2 

-cos2cpp = \ ; cos 2cpp = -5 
2cpp = 2 tt/3 cpp = tt/3 



-73- 



(4) let vp^ = tt 



-cos 2cpp 4 = -l/i = -i 
cos 2cpp 4 = 1 cpp^ = 0 
when 

vp changes by +2 TT 

cpp changes by -tt as shown below: 




I 



<s> 



E(cpp rr, k) - E(cpp,k) - -2E(k)+ E(cpp, k) - E(cpp, k) 
E($ p - n,k) - E(cpp,k) = -2E(k) 

Similarly 

F(cpp - tt, k) - F(cpp, k) = -2K(k) 

Since in this case 
\>2 ~ V]_ = 2tt 

$2 " <Pi = "TT $2 = % _ tt 

E(u 2 ) - E^) = E(cpj - tt, k) - E(cpi,k) 

E(u 2 ) - E(u 1 ) = -2E(k) 

Similarly 

u 2 - u i = -2K(k) 
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By A-43: 



l 01 ~ 



y't+i 



{+2K(k )- T ^ r 



where 



l = 



X 2 + 2 



> ^ = V 



k 2 = 2 



-L+ 1 



= + 2 + 1 



k ' 2 = 1 -* 2 =OT 




4 

(X 2 + 4) 



(£+1) = 



X 2 + 4 
2 



a ~ = A— + -j 

Jl 



t _ X 2 + 2 

n- =-^r- 



T - 4/2~ 



['(~r‘) E(k, + K(k> ] 



+ 4 



101 = 4/2 



A-53: 



6 

I 

m=l 



o 

X ^om Wi = C 01 



~ A ° Y \ 2 +2 ^ ^ XK(k) 1 

" r i/4 L Sx/PTV 00 2A?77 J 
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For n = 1: 



*12 " X 13 " *15 = 0 



By A-44: 



T _ 2{l - 1) T+2E(k) 1 

11 7ztT L— J 



21 



sm 



[ + 2K(k)] -l [Iqi] 



Since 



= 

Kl+1 J 



111 = 



' l+l 



,{+2(£+l) E(k)+ E(k) - 21 K(k) - 2tK(k)} 



*11 = 



/l+l 



{ + [( z^r) + (l+ 1} ] E(k) " 21 K(k) } 



By A-48: 



I i4 = 2 /l+l [+2 K(k) - 2E(k)] 



X 14 " 



1 + 1 



[-(*1+1) E(k) + (-t+1) K(k)] 



By A -40: 



*16 = +2tt 
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I 



m=l 



C lm I lm “ C 11 I 11 + C 14 I 14 + C 16 I 16 



X cos cp 
8/2 



tt r 



1/4 



[ I ll +I 14 + 






hi + 1 14=^== {K4 t) E(k) - <*’ K *>} 



/Z+T 



. _ x 2 + 2 . . _ X 2 . . . _ X 2 + 4 . I 2 _ (X 2 + 2) 2 

^ 2 ’ ^ 2 ’ ^ +A 2 ’ l-l 2) 2 



hi + J 14 = JJZl {+ ^^2^ E < k > - ^ K < k >} 



A -54: 



V Cl t _ A l COS$ / ■ (X 2 + 2) 2 

A m lm TTr l/4 ^ 4XA 2 + 4 
m=i 



E(k) 



-X c 



4/? 



+ 4 



t ) 



Similarly, it can be seen from the forms of C nm and Ij^ that for a com- 
plete shroud: 



I 



cos 2 $ 



2m 2m tt r 



m=l 

and in general: 
A-55: 



1/4 



[H 2 (X,k)] 



6 

I 

m=l 



Y C I 

L nm x nm 



(_l)n+l ^ cos n $ 
n 



tt r 



1/4 



[H n (X,k)] 
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where H n (X, k) is some function of X and k to be evaluated. Substituting A-53, 
A-54, and A-55 into A-26, and equating A-52 to A-26: 






A 1 cos § 



/ r (X 2 + 2) 2 ] fp . X 3 K(k) _ L tt \ 

-1/4 1 L wP7J jE(k> 77777 ^ 



n=2 



“ (-l) n A n cosn$ 



tt r 



1/4 



[H n (X,k)] 



Equating coefficients: 



A = 
o 



~ V ~ n r l/4 g 






XK(k) 

2 /x 2 + 4 



] 



A i ■ 



-V m tt r, / . a 
00 1/4 o 



H 



+ n 

4X A 2 +4 “* 4 A 2 + 4 4 



} 



A n = 0, n£ 2 



Substituting Aj into A-51: 




+Pvjri 2 rj / 4 °' 0 

4X A 2 + 4 J 4/X 2 + 4 



+ 




T(0) = 
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Recalling: 




and utilizing a lift coefficient analogous to that used in Ref. (2): 



Ct - 



T(0) 



T(0) 



(p/2)V^Lr 1/4 



p V 2 X 

r m 



2 

r l/4 



n 2 



cu 



{- 



Xi + 



r (x z + 2f 

4\/k^ + 4 



] E(k) - 



X 3 K(k) 
4 vA 2 + 4 



+ 



5 -} 



A-56: 



L 



'{[■ 



X 



T ‘ill 



X 2 + 2)' 

X A 2 + 4 



J A 2 + 4 



+ TT 



} 



This result is equivalent to that obtained by Weissinger in Ref. (1) and, as 
shown in Figure A -6 is very close to the result obtained from the more exact 
theory by both Weissinger and Morgan in Ref. (2). 



0.28 
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FIGURE A- : Lift Curve Slope as Function of Chord-Diameter Ratio. 

Taken from Prediction of the Aerodynamic Characteristics of Annular Airfoils” by W. B Morgan and E* B Caster. 

(DTMB RPT 1830) 
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